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Abstract 

The Vlasov-Einstein system describes a self-gravitating, collision- 
less gas within the framework of general relativity. We investigate the 
initial value problem in a cosmological setting with spherical, plane, 
or hyperbolic symmetry and prove that for small initial data solutions 
exist up to a spacetime singularity which is a curvature and a crushing 
Q*«\ . singularity. An important tool in the analysis is a local existence result 

with a continuation criterion saying that solutions can be extended as 
long as the momenta in the support of the phase-space distribution of 
^ 1 the matter remain bounded. 

1 Introduction 

X 

When describing the evolution of self-gravitating matter fields within the 
context of general relativity, the choice of the matter model is crucial. One 
can, for example, describe the matter as a perfect fluid, as dust, or as a 
collisionless gas. In the latter case the matter is represented by a number 
density / on phase-space, i. e., on the tangent bundle TM of the spacetime 
manifold M . The phase-space density / satisfies a continuity equation, the 
so-called Vlasov equation, which says that / is constant along the geodesies 
of the spacetime metric. Taking the energy- momentum tensor T a/3 generated 
by / as the source term in Einstein's field equations, one obtains the Vlasov- 
Einstein system for a self-gravitating, collisionless gas: 



p a GW-r^VcV/ = o, 



i 



G a > 3 = 8irT af3 , 



m 

where denote the Christoffel symbols of the spacetime metric g a p, \g\ 
denotes its determinant, G a/3 the Einstein tensor, x a are coordinates on M, 
p a the corresponding coordinates on the tangent space, greek indices always 
run from to 3, and 



m - 



1/2 



is the rest mass of a particle at the corresponding phase-space point. 

In []|] this system was investigated in the asymptotically flat, spherically 
symmetric case, i. e., with a metric of the form 

ds 2 = _ e 2„ dt 2 + e 2X dr 2 + r 2 {d0 2 + ^2 ^2 ) 

where (t,r,9,ip) are the usual Schwarzschild coordinates and A and /i de- 
pend only on t and r and vanish at r = oo. The main result was that small 
initial data lead to global, geodesically complete solutions, a result which 
has no analogue for perfect fluids or dust. Indeed, Christodoulou has shown 
that in the gravitational collapse of a dust cloud naked singularities can 
develop even for small data [Q]. Therefore, the Vlasov model seems to be 
particularly suited to describe the behaviour of matter in general relativity. 
This is further substantiated by the fact that for the Vlasov-Poisson system, 
which is the Newtonian analogue of the Vlasov-Einstein system, there are 
global existence results both in the case of an isolated system, which corre- 
sponds to the asymptotically flat case in the relativistic problem, and in the 
cosmological case, cf. [||, 12] and ||. 



In the present paper we investigate the Vlasov-Einstein system in a cos- 
mological setting. In order to simplify the problem we assume that the 
system has a high degree of symmetry and take the metric to be of the form 

ds 2 = -e 2 > x dt 2 + e 2X dr 2 + t 2 (d9 2 + sm 2 9d<p 2 ), (1.1) 

where 

{sin# for e = l, 
1 for e = 0, 
sinh^ for e = — 1, 

i>0 denotes a timelike coordinate, rS [0,1], and the functions A and fx 
depend only on t and r and are periodic in r. The angular coordinates 9 
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and ip parametrize the surfaces of constant t and r, which are the orbits of the 
symmetry action and which are spheres in the case of spherical symmetry 
e = l, tori in the case of plane symmetry e = 0, and hyperbolic planes in 
the case of hyperbolic symmetry e = — 1. They range in the domains [0,7r] x 
[0,2-71-], [0,27r] x [0,27r], or [0,oo[x [0,27r] respectively. It should be pointed out 
that the coordinates (t,r,9,cp) will in general not cover the whole spacetime 
manifold, but they do cover a neighborhood of the singularity at t = which 
will allow us to investigate the nature of this singularity. One way to think 
of the above metric is to consider the Schwarzschild metric 

2M\ 2 / 2M\ 2 2/j/i2 i • 2/ij,„2\ 



ds z = - I 1 - dt z +11- dr z + r z (d6 z + sin^ 6d^). 

If one passes through the event horizon at r = 2M then the (0,0)- and (1, 1)- 
components of the metric change sign so that the Schwarzschild radius r 
becomes the timelike coordinate. If one interchanges the notation for t and 
r and compactifies the hypersurfaces of constant t by making the components 



of the metric periodic in r one obtains a metric of the type (1J) with e = 1. 

As in || we restrict ourselves to the case where all particles have the same 
rest mass, normalized to 1, and move forward in time, i. e., / is supported 
on the submanifold 



PM:={g aP p a p f s = -l, p°>0} 



of the tangent bundle, which is invariant under the geodesic flow. Due to 
the symmetry the distribution function / can be written as a function of 
t, r, w:=e x p 1 , and F :=t 4 (p 2 ) 2 + t 4 sm^6(p 3 ) 2 . After calculating the Vlasov 
equation in these variables and the non-trivial components of the Einstein 
tensor and the energy-momentum tensor and denoting by and ' the 
derivatives of the metric components with respect to t or r respectively, 
the complete Vlasov-Einstein system reads as follows: 

dtf+ e ^L= d r f - (\w + ef-^'yJl + wt+F/tA d w f = 0, (1.2) 



er 2 ^(2t\ + l) + e = 8irt 2 p, (1.3) 
e~ 2tl (2tfi-l)-e = 8irt 2 p, (1.4) 



' = -4^ +A j, (1.5) 
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e- 2X (// + - A')) - e" 2 " (X + (A+ 1/t) (A - £)) = 8vrg, (1.6) 

where 



p(t,r):=T$(t,r) = -z / ^/l+w^ + F/t 2 f(t,r,w,F)dFdw, (1.7) 



the (3,3)-component of the field equations which is also non-trivial coincides 
with the (2,2)-component due to the symmetry. Note that on PM we can 
express p° by the other coordinates and in the above new variables obtain 



Note also that F is a conserved quantity of the geodesic flow, the modulus 
of the angular momentum of particles, and thus there is no F-derivative in 
the Vlasov equation. Furthermore, the latter equation does not depend on 



We are going to study the initial value problem corresponding to this 
system and prescribe initial data at time t=l, 



Our main result is that the solutions to this initial value problem exist on the 
time interval ]0,1] provided the data satisfy a certain smallness assumption, 
and we prove that the singularity at t = is not just a coordinate singularity, 
but a "real" spacetime singularity, a curvature and a crushing singularity. 

The paper proceeds as follows: In the next section we extract a certain 
subsystem from the full Vlasov-Einstein system (|1.2j)-( |i~6| ) and show that 
this subsystem is equivalent to the full system. In Section 3 we prove a local- 
in-time existence and uniqueness result for classical solutions, together with 
a continuation criterion which says that when going backward in time, i. e., 
towards the singularity a solution can be extended as long as the support 
of f remains bounded with respect to w. The main difficulty here is to 
show that a solution cannot break down due to a blow-up of a derivative of 





e. 



f(l,r,w,F) = f(r,w,F), A(l,r 



o 

)=A(r), ^(l,r)=/l(r). 
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/ (or one of its moments), that is to say, there is no formation of shocks. 
This continuation criterion is used in Section 4 to prove that solutions exist 
on ]0,1] for sufficiently small initial data, the support of / with respect to 
w is shown to decay like t c for some c>0 as t— >0. The structure of the 
singularity at t = is analyzed in Section 5. In the last section we briefly 
investigate the behaviour of the solutions for t>l. In order to extend a 
solution forward in time one needs to bound the support of / with respect 
to w and the metric component e 2>l . A bound on the former quantity can 
be established regardless of the size of the initial data, but for e = 1 it can 
be shown that e 2/i blows up in finite coordinate time. 

As to the physical relevance of the situation studied here our point of 
view is that some results in general relativity are intended to describe con- 
crete, real-world phenomena while others are intended to elucidate general 
features of the theory like for example the structure of possible singularities. 
The present paper belongs to the second category. 

An approach which may generalize to other situations more easily than 
the present one but gives less information on the structure of the singularity 
and on questions of existence of solutions is taken in 11], where the system 
with spherical and plane symmetry is analyzed using constant mean cur- 
vature slicing. As is shown in |l(J homogeneous solutions, i. e., solutions 
which are independent of r have a curvature singularity provided / is not 
identically zero. 

To conclude this introduction we mention some further results on the 
asymptotically flat, spherically symmetric Vlasov-Einstein system: In || it 
is shown that solutions of this system converge to solutions of the Vlasov- 
Poisson system in the Newtonian limit. In || it is shown that if a singularity 
forms the first one has to form at the centre of symmetry. The existence of 
static solutions is established in []|, [?]]. 

Acknowledgements: I would like to thank A. D. Rendall for helpful 
discussions and comments. The present investigation was completed during 
a stay at the Erwin Schrodinger International Institute for Mathematical 
Physics in Vienna. I would like to thank the Institute and in particular 
P. Aichelburg and R. Beig for the kind invitation. 



2 Equivalent subsystems 

Let us first make precise the regularity properties which we require of a 
solution: 
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Definition 2.1 Let /ClR be an interval. 

(a) /€C 1 (/x]R 2 x]R+) is regular, if f(t,r + l,w,F)= f(t,r,w,F) for 
(t,r,w,F) G / x H 2 x Hq", />0 , and supp/(t,r,-,-) is compact, uni- 
formly in r and locally uniformly in t. 

(b) p (or p, j, q) gC 1 (/x]R) is regular, if p(t,r + 1) = p(t,r) for (t,r) G 
IxE 

(c) AeC^JxE) is neuter, i/AeC^IxM) and A(t,r + 1) = \(t,r) for 
(t,r)eIxM. 

(d) ^.gC^/xIR) is regular, if p! G C 1 (7 x H) and p(t,r + l) = p(t,r) for 
(t,r)eJxIR. 

We identify such functions with their restrictions to the interval [0,1] with 
respect to r. The fact that regularity means different things for different 
objects will cause no ambiguities. 

Let (f,X,p) be a regular solution of the subsystem (|L"2j), (|Q|), (|1.4| ) on 
an interval / with 1 G /. We want to show that ([O]) and hold as well. 
Integrating (|1.4| ) we obtain 




(2.1) 



and 




From ([l]^) and integration by parts it follows that 



a 



i 

Jo 



oo 



W 



,2 



p'(s,r)s 2 ds = 7r 



o y/l + w 2 + F/s 



=d r f(s,r,w,F)dFdwds 



— oo 




Adding ((T|) and yields 



A + / u = 4vrte 2/ "( / 9+p) 



(2.2) 
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and if we assume that the constraint equation Q1.5Q holds at time t = 1 then 
these identities imply 

te~ 2 ^ (// + 4yr/;e A+M j) = -4vr j\p+P) (y«'+47rse A+ ^ s 2 ds 

so that 

/ u / + 4vrte A+/i j = 

on /, i. e., (|l.5|) holds for all iE/. The latter equation can be differentiated 
with respect to r to yield 

/ u" = (A' + / u> , -47rte A+ ^j'. 

From (IO) we obtain by integration by parts the identity 



f(t,r) = 



oo roo 



oo JO 



l + w 2 + F/t 2 d t f 
+ ( XwJl + w 2 + F/t 2 + e' 1 - x f i\l + w 2 + F/t 2 ) ) d w f 



dFdw 



/OO /*00 y 
/ a/1 + w 2 + F/t 2 d t fdFdw - e x -v\( p +p) - 2p! j. 
-ooJO V 



Equation (1.3) can be rewritten in the form 



\ = ^te 2 »p-^—. (2.3) 



Since 



2p(t,r) vr [°° f°° F/t 3 
P(t,r) = / / —======fdFdw 



t t 2 „ '-co Jo \f\ + W 2 + F/t 2 ' 

y ^ji+w 2 +F/t 2 d t fdFdw 



2p(t,r) 2q(t,r) , tt 



-oo JO 

differentiating ( |2.3| ) with respect to i yields 



7T Z" 00 Z" 00 / 

+ J \/l + w 2 + F/t 2 d t fdFdw, 



A = _W 1 p + 2A/i+^-87re 2 ^ 

+ ^ e 2„ /• L . o . ^,.o „ , ,., , 1+--" 



/ „ 1 y ^/i+w 2 +F/t 2 d t fdFdw , o/2 . 

Combining these identities implies the remaining field equation (|1.6| ). Thus 
we have established the following result: 
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Proposition 2.2 Let (/,A,/x) be a regular solution of (|L2|), (|L3|), ( |1.4D on 

some tame interval I C IR + 1 € I, and £/ie initial data satisfy (|L5|) 
fort = l. Then (|L5|) and (|ll| /joW /or a// t G /. 



When proving local existence it now suffices to consider the subsystem fll.z 
( |1.3| ), Jl.4| ). However, it would then become technically very unpleasant to 
control // . It is more convenient to consider an auxiliary system, which 
consists of the modified Vlasov equation 



d t f+ Vl + J™ F/f , d rf ~ (Aw + e^ x ^l + w^ + F/t^ d w f = 0, (2.4) 
together with (|0|), ([yD, and 

f L = -4iTte x+fl j. (2.5) 

Assume that we have a regular solution (/,A,/x,/2) of this system, where 
regularity for jl means that this function has the same properties as // for 
fj, regular. We want to show that \x is nothing else than jj so that by 
Proposition (/, A,/i) solves the full system. As above, 

tfj,'(t, r)e- 2 » = /j'(r)e~ 2M ° + ^ j\\s,r) s 2 ds 

and 



p' s 2 ds = - I (\ + fi)e x ~ fM js 2 ds-e x ~Vj s 2 -j fi{p+p)s 2 ds. 



s=t 



Using ( p^2; ) and (2J3) we obtain 

M 't e - 2 ^ = tf e - 2 $ - A7rt 2 e x -^ j + ATie x -^j 
so that ( |T75|) holds for all i € / if it holds for t = 1. 



Proposition 2.3 Lei (/,A,//,/i) 6e a regular solution of fl2.4|) , ( |l.3|) , (|T 

Q. T/ien (/,A,/x) sofees 

We conclude this section with a result which reflects the conservation of the 
number of particles in our system and is an immediate consequence of the 
Vlasov equation. 



8 



Proposition 2.4 Assume that f is regular and satisfies Q1.2D with regular 
coefficients A and \i. Then 

rl poo roo 

/ / / e x ^ f{t,r,w,F)dFdwdr 

JO J-ooJO 



and 

"1 ('OO /■oo 



is conserved. 

3 Local existence and continuation of solutions 

In this section we prove the following local existence and uniqueness result 
with the continuation criterion described in the introduction: 

Theorem 3.1 Let fe C^IR 2 x H+) with f(r + l,w,F) = f(r,w,F) for 
(r,w,F)eM 2 xIR+, />0, and 

wq : = sup | (r,w,F) € supp/| < oo, 
Fq := sup j.F | (r,w,F) G supp/j < oo. 

Let A^eC^H) with A(r) = A(r + l), jl(r + 1) = for rGE and 

/ 5 / (r) = -4vre A ° +/! j(r), rGlR. 

In the case of hyperbolic symmetry e = — 1 assume in addition that /5(r) < 
/orr€lR. T/ien there exists a unique, left maximal, regular solution (f,X,fi) 



of (1.2)-(|1.6|) with (/, A,/x)(l) = (/, A,/5) on a time interval ]T,1] with T€ 
sup||w| | (t,r,w,F) Gsupp/| < oo 



[o,i[. If 



then T = 0. 

Remark: To motivate the restriction on jl in the case e = — 1 consider the 
following "pseudo-Schwarzschild" solution 

ds 2 = _ fi + 2 J± ^" 1 d f + (l + ^£ W + 1 2 ( ^2 + ginh 2 ^2) 
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with M € IR which is a vacuum solution of our system. The restriction jl< 
is equivalent to M > so that this case is contained in the present investi- 
gation. For M = the spacetime is flat, and for M<0 it has a coordinate- 
singularity at t = — 2M, but the spacetime can be extended through this 
singularity, something we want to exclude from our investigation. 
Proof of Theorem 3.1: Define 

S a, 

and consider the auxiliary system ( |2.4[ ), Ql.3| ), ( |1.4|) , (|2.5|). We construct a 
sequence of iterative solutions in the following way: 

o o 

Iterative scheme: Let Ao(t,r) := A(r), fi (t,r) := fi(r), fL (t,r) := fi{r) for 
]0, 1], r € IR. If A n _i, /U n _i, /i„,-i are already defined and regular on ]0, 1] x ]R 
then let 



G n -x(t,r,w,F) :-- 



pf-ln — l -^n — 1 jii I 



\y/l + W 2 +F/t 2 ' , 

and denote by (R n , W n )(s,t,r,w,F) the solution of the characteristic system 

-^-(R,W) = G n „ 1 (s,R,W,F) 
as 

with initial data 

(R n ,W n )(t,t,r,w,F) = (r,w), (t,r,w,F) e]0,l] x IR 2 x IR+; 
note that i* 1 is constant along characteristics. Define 

f n (t,r,w,F):=f((R n ,W n )(l,t,r,w,F),F), 
that is, f n is the solution of 

g/^n — 1 Aj-j, — 1 qjj 

dtfn+ y /i + w * + F/* drfn 



\ n ^w + e^- A ^fl n ^l + w 2 + F/t 2 ) 8 w f n = 

with /n(l) = / ; an d define p n , p n , j n , q n by the integrals (P1)- (|i~Tol ) with 
/ replaced by f n . Recall that the solution of ([h^) is given by ([O]) so we 
define \i n by 

e -^( t „) = f^±l_ e _^/' p „ (S|r)s2ds . p.!) 
t t Jl 
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note that the right hand side of this equation is positive on ]0,1] x IR. It is 
at this point that we need our additional assumption on the initial data in 
the case e = — 1. Finally define 

X n (t,r):=47Tte 2 ^ Pn (t,r) (3.2) 

cf. Q, 

X n (t,r):=X(r)+ \ n (s,r)ds, 



fl n {t,r) := -47rte^ +x -j n (t,r). (3.3) 

These iterates are regular on the time interval ]0,1], in particular, since 
A„_i, fJL n -i, A n _i, jjL n -i are continuous on ]0,1] x IR and periodic in r, these 
functions are bounded on compact subintervals of ]0, 1], uniformly in r, and 
since G n -\ is linearly bounded with respect to w the characteristics R n , W n 
exist on the time interval ]0,1]. 

The proof of Theorem |3.1| now consists in showing in a number of steps 
that the iterates constructed above converge in a sufficiently strong sense. 
Step 1: As a first step we establish a uniform bound on the momenta in the 
support of the distribution functions f n , more precisely we want to bound 
the quantities 

P n (t) :=sup{|u;| I (r,w,F) Gsupp/ n (t)| 
uniformly in n. On supp/ n (t) we have 

l + W 2 +F/t 2 < v /l + P n (t)2 + jFo / t 2<i±^0 (1 + p n(t))) 



and thus 



and 



(1 + Fn) 2 ° 



\\Pn(t)\\, ||jn(*)||<c§||/||P n (*) 2 . 



t z 

Throughout the paper || • || denotes the L°°-norm on the function space in 

o 

question; we have used the fact that ||/n(*)ll = ll/ll for and ie]0,l]. 

The numerical constant c may change from line to line and does not depend 
on n or t or on the initial data. In view of the continuation criterion it 
is important to keep track of any dependence on the latter. From (|3.1D it 
follows that 

e -2/i n (t,r) > £1 
~ t 
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where 



ci = ci(/2) := 



infe 2 ^ for e = lore = 0, 
mle~ 2 ^ — 1 for e = — 1. 



By (^2), (3J3), and the above estimates on p n and j n we get 



An(s,r) <47rse 2 ^|j n (s,r)|<c- 



ci 



l^n(s)' 



and 



A n (s,r) < Airse 2 ^p n {s,r) + 



l + e 2 ^ 
2s~ 



Thus 



where 



<< (1 +^ii/i|(I±^ + I^i 

ci s 2s 



W n+1 (s) <-(l + P n (s)) 2 (l + \W n+1 (s)\), 



C2 = c 2 (/,Fo,/l):= C (l + l/ C i)(l + Fo) 2 (l + ll/H). 
This implies that 

Pn+l(t)<W + C 2 -(l + P n ( S )) 2 (l + P n+ i(s))ds. 

Jt s 

Let z\ be the left maximal solution of the equation 

zi(t) = w + c 2 / -(l + zi(s)) 3 (is, 
Jt s 

which exists on some interval with T\ £ [0,1[. By induction 

Pn(t)< Zl {t), t€]Ti,l], nGM, 

and all the quantities which were estimated against P n in the above argument 
are bounded by certain powers of z\ on 

Step 2: Here we establish bounds on certain derivatives of the iterates. In 
particular we need a uniform bound on the Lipschitz-constant of the right 
hand side G n of the characteristic system in order to prove convergence in 
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the next step. Differentiating ( |3.1| ) and ( |3.2| ) with respect to r one obtains 
the identities 



p 2 Vn / o ft \ 

&(t,r) = — (y?(r)e- 2 " + 47r ^ p' n (s,r) s 2 ds) , 

\' n (t,r) =e 2 ^ ^87r^(t,rK(t,r)+4vr^(t,r)-^;(t,r)), 

A^^r) = A / (r)+/ t A;( S ,r)d S . 



In the following C\ denotes a continuous function on ]Ti,l] which depends 
only on Z\. By Step 1, 

ll/4(f)ll> IlKWII, <Cl(*)ll^-/«(*)II- 

Define 

L> n (t):=sup{||a r / n (s)|||i<s<l}. 

Then the above estimates and the formulas for the derivatives of the metric 
components show that 

K(t)||, ||A^(t)||<C 1 (i)( C3 + L' n (t)), 
where c 3 := ||e _2/ V|| + ||A'|| + 1. From it follows that 

e^- x "jl n = -ATTte 2 ^j n , 

and 

(e^- A »/i n ) ' (t , r) < Ci (t) (c 3 + D n (i) ) • 

We are now in the position to estimate the derivatives of G n with respect 
to r and w: 



d r G n (t,r,w,F) = (( Mn -A n )'e^- A " 



W 



sJl + w 2 + F/t'^ 
-(e^- K fi n y^l + w 2 + F/t 2 -X' nU ?j . 



d w G n (t,r,w,F) = (e^~ x - 



^l + w 2 + F/t^ , 



cP'TL Sin . , 



y/l + W 2 + F/t 2 J 
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and thus 

\d r G n (t,r,w,F)\ < Ci(t)(c3 + D„(t)), 
\d w G n (t,r,w,F)\ < Ci(t) 

for f€]Ti,l], r€lR, Fg[0,F ], and |w|<zi(t). Differentiating the charac- 
teristic system we obtain 

-^-d r (R n+1 ,W n+1 )(s,t,r,w,F) < 

Ci(s)(c 3 + D n (s)) \d r (R n+1 ,W n+ i)(s,t,r,w,F)\, 
and thus for (r,w,F) € supp/ n+ i(t) Usupp/ n (t) 

|9 r (i? n+ i,^ n+ i)(l,t,r,u;,F)|<exp ^ C\(s){c 3 + D n (s)) ds^ 
By definition of D n this implies that 

D n+1 {t) < \\d {r:W) f\\ exp Q^CiW (c 3 + A»(a)) ds) . 
Let ^2 be the left maximal solution of 

25 (*) = l|d(ryu,)/l|exp Ci(fi) (c 3 + Z 2 ( S )) ^ 

which exists on an interval ]?2,1] c]Ti,l]. Then by induction, 

D n (t)<z 2 (t), t€]T 2 ,l], neJN, 

and all the quantities estimated against D n above can be bounded in terms 
of z 2 on ]T 2 ,1]. 

Step 3: Let [<5, 1] C]T 2 , 1] be an arbitrary compact subset on which the es- 
timates of Steps 1 and 2 hold. We will show that on such an interval the 
iterates converge uniformly. Define 

a n (*):=sup{||/ n+ i(r)-/ n (r)|||r€[*,l]}, 

and let C denote a constant which may depend on the functions z\ and z 2 
introduced in the previous two steps. Then 

\\pn+l(t)-p n {t)\\, \\Pn+l(t)-p n (t)\\, \\j n+1 (t) - j n (t) \\ < Ca n (t) , 
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and thus 



Therefore, 
By Step 2 



|A n +i(i) — A n (£)||, ||A n+ i(t) — A n ( 
WfJ-n+i(t) - fJ. n (t)\\, \\p, n +i(t) - p, n (t)\\ <Ca n (t). 

\G n+ i - G n \(s,r,w,F) < Ca n -i(s). 

<C 



d( r ,w)Gn(s,r,w,F) 

for all sE[S, 1], nGlN, and (r,w,F) with \w\ < zi(s). For characteristics 

o 

which start in supp/ this implies 



^-(R,W) n+1 --^-(R,W) n 



(s,t,r,w,F)< 
c\(R,W) n+1 -{R,W) n \(s,t,r,w,F) + Ca n - 1 (s), 
and by Gronwall's inequality 

(R,W) n+1 - (R,W) n \(l,t,r,w,F) < C j\ n ^{s)ds. 
If we recall how /„ was defined in terms of the characteristics this implies 

®n(t)<cj^ a n _i(s)ds, n>l. 
By induction we obtain 

, , „C n (l-t) n C n+1 



for n G IN and t G [5, 1] . This implies that f n and all the other quantities 
whose differences were estimated in terms of a n converge on [6, 1] , uniformly 
with respect to all their arguments. These quantities therefore have con- 
tinuous limits, but the established convergence is not yet strong enough to 
conclude the differentiability of, say, / :=lim n -, 00 f n . In order to achieve the 
latter we need the following lemma: 

Lemma 3.2 Let (A,/x,/i) be regular on some interval /ClR + , and let 
(R,W)(-,t,r,w,F) be the solution of 



n-\ 



w 



r = 



y/T+uF+Fjs 2 ' 



w - 



-\(s,r)w-e fM ~ x fi(s,r)^l + w 2 + F/s 2 
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with (R,W)(t,t,r,w,F) = (r,w) for (t,r,w,F) € I x IR 2 x Mj" . De/me 

f( a ) := e^-^^^s.^r,^,^), 
:= 9H^(s,t,r,w,F) 



+ a/I + ^ + F/sV^A dR(s,t,r,w,F) 

/or (?G {9 r ,9tu}. T/ien these quantities satisfy the following system of differ- 
ential equations 

i(s) = a 1 (s,R(s),W(s),F)as) + a 2 (s,R(s),W(s),F)i 1 (s), 

17(a) = (a 3 + a 5 )( S ,R(s),W(s))( S ),F)Z(s) + a 4 (s,R( S ),W(s),F) V (s), 



where 



w 2 



ai{s,r;w,F) := — — 2 A-/t, 



a 2 (s,r,u;,F) 



1 + F/s 



2 



(l + w 2 + F/s 2 ) 3 / 2! 



a 4 {s,r,w,F) := — — =— M* fi+ ==A , 



a 5 (s,r,w,F) := -yjl + w 2 + F/s 2 e 2 ^ 

(e- 2A (p! + £(//- A')) - e" 2 " (X + (A + (A - £))) . 

In particular, ifp, = fj,' and (A,/i) solves the field equations ( p..3| )-( [L^ ) t/ien 



a,5(s,r,w,F) = — y l + u> 2 + F/s 2 e 2 ^ r) 8Trq(s,r) . 

The proof is only a lengthy calculation and therefore omitted. The lemma 
will be used twice in the further argument: In the next step it will be used 
to prove that also certain derivatives of the converging sequences obtained 
in the previous step converge, thus obtaining a regular, local solution. Then 
it will be used in Step 6 to show that control on the support of the solution 
with respect to w suffices to extend the solution to the interval ]0,1]. 
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Step 4 : Fi x ^£]^2)1] and U>0, and consider the system derived in 
Lemma 3^2 with (A n ,/x n ,/i n ) instead of (A,/i,/i), and call the correspond- 
ing coefficients a U; i, i = l,...,5. By Steps 1 and 2 we have the estimates 

\a n> i(t,r,w,F)\+ d^^a n< i(t,r,w,F) <C (3.4) 

for n £ IN, i = 1,. . . ,4, < i 7 < i*b, < £/, and t € [(5, 1]. The only new terms 
to estimate here are fi n and fi' n , but from ( |3.1| ) we obtain 

^ n = Airte^pn + , 

and 

An = 2/4 (An " - ) +47Tte 2 ^K, 

so both of these terms are bounded by Steps 1 and 2. The convergence 
etablished in Step 3 shows that 

a n> i(t,r,w,F)-a mt i(t,r,w,F)-*Q, ra,m^oo, i = l,...,4, 

uniformly on [5,1] x IR x [ — C7, t/] x [0,-Fo]- The crucial term in the present 
argument is a n ^, more precisely the expression 

9n := e~ 2Xn (fi' n + jln{^n ~ K)) ~ ( A n + (A n + l/t){X n - fl n )) . 

If the iterates solved the field equation ( |1.6| ) then this term would equal 
8irq n and would also converge. The idea how to treat a n ^ is to show that 
g n — 8irq n — > for n — > oo and then use the fact that g n converges and has 
uniformly bounded r-derivative. Now 

An = (/4 + A„)An " 47Tte^ +A " j; 

and 

A n = 2 A n A„ + £p + 47rte 2 "» ^n + ^) + ■ 
From the definition of p n we obtain 
. 2 2 
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where we used the Vlasov equation to express dtf n and integrated by parts; 
note that the coefficients in that equation have index n — 1 . Inserting all 
this into the expression for g n yields, after cancelling a number of terms, 



pn— 1 



g n = 2e ""'fin \ji; n -e 

+4:TTtj' n (e fln - 1 ~ Xn - 1 -e 
+e~ 2 ^ (A„ + An)(A n -i - A n ) + 87rg n . 



— Ar 



By Steps 1, 2, and 3 it remains to show that \J n — /2n-i^0 for n— >oo in 
order to conclude that g n ^ 8irq n . To see the former differentiate ( |3.1| ) with 
respect to r to obtain 



Mr, 



Differentiating the defining integral of p n , using the Vlasov equation for f n 
to express d r f n , and integrating by parts with respect to w and s results in 
the relation 



Mn : 



-p!e- 2 t + 47re x -^) +te ^-^„-i-A n +A n _ 1 



+ - 



Mr, 



1 (A n -l + An-l)An ~~ (\i + An)An-l 



(is, 



Mr 



• /i for n — ► oo, in 



and since the initial data satisfy the constraint ([[ 

particular, /i^ — Mn-i — ► for n — > oo. 

In Step 3 we have shown that among other quantities the characteristics 
(Rn, W n )(l,t,r,w,F) converge. Now for any e>0 there exists iVeM such 
that for n,m > N, s G [6, 1] , r 6 IR, \w\ < U, and F 6 [0,Fq] we have 



a n , 5 (s,r,u;,F) - (-8vr)e 2 ^ yjl + w 2 + F/s 2 q n (s,r) 



<£ 



and 



|d r g n (s,r)| <C, jg n (s,r) -g m (s,r)[ <e, 
which together with the estimates ( |3.4[ ) implies that 

£n - Cm (s) + | fjn ~ Vm | (s) < C £ + C |£ n - £ m | (s) + C | r/ n - 7/ m | (s) . 

This implies the convergence of d( r ^(R n ,W n )(l,t,r,w,F); note that the 
transformation from d(R,W) to in Lemma |3.2| is invertible, and the 
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coefficients in the transformation are convergent in the present situation. 
Thus the limiting characteristic (R,W)(l,t,r,w,F) and therefore also / are 
continuously differentiable with respect to r and w. This in turn implies 
that all the moments calculated from / are differentiable with respect to 
r, the coefficients in the limiting characteristic system are continuously dif- 
ferentiable with respect to r,w, and F, and thus (R, W)(l,t,r,w,F) is also 
differentiable with respect to F and t, and the regularity of the limit (/, A,/i) 
is established. 

Step 5: The estimates on the difference of two consecutive iterates derived 
in Step 3 can also be used on the difference of two solutions / and g with 
the same initial data. This results in the estimate 

sup{||/(r)-^(r)|||r€[t J l]}<Cjr 1 sup{||/(r)-^(r)|||rG[ Sj l]}d S 

on any compact subinterval of ]0, 1] on which both solutions exist, and thus 
f = g there. 



Step 6: To conclude the proof of Theorem 3.1 it remains to establish the con- 
tinuation criterion. Let (f,X,/j,,p,) be a left maximal solution of the auxiliary 
system (fOj) , (|L3|) , (jT|), with existence interval ]T,1]. By Proposi- 



tion (/,A,/x) solves (|1.2|)-( |L6| ). Now assume that 

P* :=sup{H \(r,w,F) £supp/(i), t€]T,l]} < cx>. 

We want to show that T = so let us assume that T > and take to G 
]T, 1[. We will show that the system has a solution with initial data 
(f (to) , \(to) , /j,(to)) prescribed at t = to which exists on an interval [to — 5, to] 
with 5>0 independent of to- By moving to close enough to T this would 
extend our initial solution beyond T, a contradiction to the initial solution 
being left maximal. 

Steps 1-5 have shown that such a solution exists at least on the left 
maximal existence interval of the solutions (^1,^2) of 

to 1 

s< ds, 



Zl (t) = Wo + c 2 [ -(l + zx(s)) 
Jt s 



Z2{t) = \\d( r)V) )f(t )\\exp(J Ci(s)(c 3 + z 2 (s))ds 

where 

W :=sup||w| I (r,w,F) esupp/(i )|, 
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ci = ci(fi(t )):-- 



infe 2 M*o) for e = lore = 0, 
i n f e -M*o)_i for e = -l, 



C2 = c 2 (/(to),i ? o,/u(io)):=c(l + l/c 1 )(l + Fo) 2 (l + ||/(to)||), 
C 3 = ||e- 2 ^') / u'(to)|| + ||A'(t )|| + l, 

and the function d depends on z x . Now W <P*, \\f(t )\\ = ||/||, F 
is unchanged since -F is constant along characteristics, and ( |2.lD shows 
that ci(/u(to)) > ci(/2). Thus there exists a constant C2>0 such that 
c 2 (/(t ),i ? o,^(io))/s<c^ for t £]T,l] and sg[T/2,1]. Let z\ denote the 
left maximal solution of 

zl(t)=P* + c* f t \l + z* 1 (s)fds. 



Next observe that the coefficients a± , . . . , 05 in Lemma p.2| are uniformly 
bounded ]T, 1] along characteristics in supp/ if we let (i = fi' and use the 
field equation ( |1,6| ). The lemma then shows that 

zr :=Bup{||a (ritB) /(t)|||r<t<i} <oo. 



From 



A'(t,r) = e 2 ^ ^87rt/i' (t, r)p(t,r) + 47rtp'(t,r) - ^//(t,r)) , 
A'(t,r) = A'(r) + f \'(s,r)ds 



we obtain a uniform bound C3(/i(to), A(to)) < C3. Let z| be the left maximal 
solution of 

z\ (t) = D* exp Qf *" C{ (s) (4 + z*(s))d s y 

where depends on z\ in the same way as C\ depends on Z\. Clearly, 
z\ and z\ exist on an interval [to - 5, to] with <5>0 independent of to- If we 
choose 5<T/2 then z\<z* and z 2 < z 2 by construction, in particular, zi 
and z 2 exist on [to — 5, to] , and the proof of the continuation criterion and of 



Theorem 3.1 is complete. □ 
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4 Existence up to t = 

In this section we show that the solutions obtained in the previous section 
exist on the interval ]0, 1] provided the initial data are sufficiently small. 



Theorem 4.1 Let (/, A,/t) be initial data as in Theorem |3.1| , and assume 

that 

c := I _ io ttVjFo ^/l + w 2 + F \\e 2 ^ \ \ \\f\\ > 

in case e = or e = 1, and 

2° 

c:=i (l- He 2 !) - 107r 2 Wo F ov /l + W g + Fo ^^JL ||/|| > 



^n case 



9 ll 

e= — 1 — note that \\e ^\\ < 1 in this case by the assumption in Theo- 
rem 3.1 . T/ien i/ie corresponding solution exists on the interval ]0,1], and 

\w\<w t c , (r,w,F) Gsupp/(i), ie]0,l]. 

o o 

Proof : Let (/,/!, A) be initial data satisfying the smallness assumption, and 
define 

P(t) :=sup{|u;| | (r,w,F) Gsupp f{t)}, te]T,l]. 

The characteristic system for (|l.2| ) and the field equations (|l.3| ) and (|l.5| ) 
imply that 



u, = -Aw-e^VVl + wHF/i 2 

l + ee 2 » 



Ante 2 ^ (^j^l + w 2 + F/t 2 -pu?j +■ 



2i 



-u>. 



Assume that P(t)<wo for some f€]T,l], which is true at least for i = l. 
Then 

0<p(t,r)<^ / Jl + w 2 + F/t 2 f(t,r,w,F)dFdw 
t 1 J- Wo Jo v 

< 2^0^0/1+^+^ 11/11 1~ 3 , 

and 

t jO JO 2 



t 2 J-P(t)Jo ' w ' ' ' " 2 
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Next we have the estimate 



e -2M*,r) = f ±l_ e _8vr / S 2 p( S ,r)d S >- — -€ 



so that 



where 



e -2 At (t,r) > £1 

Jinfe -2 ^ for e = 0ore = l, 
Cl ' |infe~ 2/? -l for e = — 1. 



Thus e 2 ^*' r ) < cjH, and 



1 + ee 2 ^ 1_C2 
2t ~2t~' t 



for e = or e = 1, 



1 + 1 / 1 \ 1 q 1 c 2 



2t ~ 2t V 1 + cit- 1 / ~ 2 1 + ci t 't 
for e = — 1. Assume that w(t)>0. Then 

1 SiA w (0 



w(t)> [c 2 -87r 2 woF oy Jl + w% + Fo 



-2ir 2 w F ^l + w 2 + F ±- \\f\\ (4.1) 



whereas if w(t) < we obtain the estimate 



w(t) < [c2-8tt 2 w F J1 + w^ + F 



1 ° A w(t) 



Cl / t 



1 „ h Pit) 



+2ir 2 w F Jl + w 2 + F -\\f\\^. (4.2) 



If we let t = 1 in ( t4,l|) and (|4.2| ) our smallness condition on the initial 
data implies that there exists a small constant <5>0 such that w(l)>0 if 
wo/(l + 5) < w(l) < wo, and w(l) < if — u)o < ^(1) < — wo/(l + ^)- This im- 
plies that P(i) < i«o ° n some interval ]to,l[ which we choose maximal with 
this property. On the interval ]io,l] the estimates fl4.l|) and (|4.2|) hold for 
any characteristic which runs in supp/ and for which w(t) > or w(t) < 
respectively. 
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Let iE]io)l] be such that w(t)>0 for a characteristic in supp/, and 
choose t\>t maximal with w{s) >0 for s£ [t,t\[. Then fl4.1|) holds on [t,t\[ 
which by Gronwall's inequality implies that 



w(t) < exp ^c 3 y — 
w 



w{ti)-Ci / exp( -c 3 ^ — 



(ti) + c 4 ti C3 J\- l - C3 P{s)ds 



dr \ P(s) 



ds 



where 



c 3 :=c 2 -8vr^oi ? oV 1 + u, o+ i? o — ll/ll, 

Cl 



c 4 := 2vr^ i ? oV 1 + u; o+ i? o-| 

^1 



If = 1 then 



u>(i)<i C3 ^wo + c 4 j\- 1 - C3 P{s)ds Ss j . 
If £i < 1 then w(ti) = 0, and 



W(t) < (t/tl) C3 C 4 ^ 3 



-1-C3 



P(s)ds 



< t C3 ^w + c 4 j\~ 1 - C3 P(s)ds 



Consider now ie]io,l] such that w(t)<0, and choose ii>i maximal with 
w(s) <0 for s€ [t,ti [. Repeating the above argument, but now using 
instead of fl4,l|) , yields the estimate 



w(t)>(t/h) ca ^(ti)-c 4 ti C3 y s~ 1 - C3 P(s)dsj, 
and distinguishing the cases ii = 1 and t± < 1 as above implies that 

<* C3 (wo + c 4 j\- l - C[i P(s)ds^j . 



Therefore, 



P{t)<t cz [wQ + c 4 j\- l - C3 P{s)ds 
for all te]to,l]- Applying Gronwall's inequality again yields the estimate 

P(t)<w t c *- C \ te]to,l], 
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since z(t) :=wot C3 ° 4 is the solution of the integral equation 
z(t) = t° 3 (wQ + aj s~ l - cz z{s)ds\ 

which is equivalent to the initial value problem 

z(t) 

Z(t) = (c 3 -C4) — , Z(1)=W . 

The estimate on P(t) implies in particular that to = T, i. e., it holds on the 
whole existence interval of the solution, which by Theorem [O] implies T = 0, 
and the proof is complete. □ 



5 The nature of the singularity at t = 

In this section we investigate the behaviour of solutions as t— >0. Theo- 
rem 4.1 shows that there are solutions which exist on ]0, 1] so that we do not 
investigate the empty set. However, the results which follow do for the most 
part not depend on the smallness assumption which we used in the previous 
section. First we show that solutions which exist on ]0,1] have a curvature 
singularity at t = 0, more precisely: 

Theorem 5.1 Let (/,A,/z) be a regular solution of ( |l.2D -( |l.6D on the inter- 
val ]0, 1]. Then 

(R a p l5 R a ^ s ){t,r)>Q{inie~ 2 ^ + e) t~ 6 , t€]0,l], reH, 

where R a ^ S denotes the Riemann curvature tensor corresponding to the 
metric given by A and \x. 

The curvature scalar considered here is sometimes called Kretschmann 
scalar. 

Proof : It can be shown that 

R a ^sR al3yS = 4( e - 2A ( / x" + / u / (/i'-A / ))-e-^(A + A(A-/ l ))) 2 



_8 



K x + K 2 + K z . 
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One way to see this is to use a computer algebra system like Maple, another 
is to exploit the symmetry of the metric and split the summation 



AbCd 



Ruf3~i5R al3 ~' 6 = RabcdR abcd + RaBCdR ABCD + ^RAbCdR 

where lower case Latin indices take the values and 1, and upper case Latin 
indices take the values 2 and 3. Now 

Rabcd = (e- 2X (^ + fi\fi'-X f ))-e- 2 ^X + X(X-fi)))(g ac g bd -g bc g ad ), 
Rabcd = t~ 2 (V 2 ^ + e) {gAcQBD - 9bc9ad) , 

RAbCd = t~ l 9AC^bd^ 

which can be seen after a lengthy calculation. Inserting these expressions 
into the above summations yields the formula for the Kretschmann scalar. 
The first term K\ is clearly nonnegative and can be dropped. In order to 
estimate K 2 we insert the expressions 

e + e" 2 ^ 



e-^A = Airtp 
e~ 2fl p = Antp + 



2t ' 

e + e- 2 ^ 



2t 



into the formula for K 2 and obtain 
+2 

—K 2 = 16ir 2 t 2 (p 2 +p 2 -2j 2 )- 47rt(p-p) 



t 2 ~ 2 2, 2 2 2^ „ , ,e + e"^ (e + e" 2 ^) 2 

t 2t 2 ' 



Now 

< p(t,r)'' 2 p(t,r)'< 2 
by the Cauchy-Schwarz inequality. Therefore 

p 2 +p 2 - 2j 2 > p 2 +p 2 - 2pp = (p-p) 2 

and 

t 2 , e + e' 2 ^ (e + e- 2 ^) 2 

-K 2 > (4 7r t(p-p)) 2 -4irt(p-p) e -±^+ [€ + e 2t2 > 

/ , . e + e- 2 ^\ 2 l(e + e- 2 ^) 2 
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Thus 

2 



f e + e -2Ai(t,r)V 

(R aM5 R a W) (t,r) > 6 ^ ¥ L 



Recalling ( J27l| ) we obtain 

e M + e = / s z p(s,r)ds> , 

t £ A ~ t 

and inserting this into the previous estimate completes the proof. □ 

Next we show that the singularity at t = is also a crushing singularity in 
the sense that the mean curvature of the surfaces of constant t blows up as 
t->0. 

Theorem 5.2 Let (/, A,/i) be a solution of ( |1.2|) -( |F6| ) on ]0, 1] and define 
c:=infe~ M for e = or e = 1, andc:= |inf(e _2M — l) 1 / 2 /or e = — 1. Lei fc(i,r) 
denote the mean curvature of the surfaces of constant t. Then 

k{t,r)<-ct' 3/2 . 

Proof : For a metric of the form 

ds 2 = —ip 2 (t,x)dt 2 + gijdx l dxi 

where i,j run from 1 to 3 the second fundamental form is given by 

kij = -(^y 1 9ij, 

and its trace k(t,x) = k\{t,x) is the mean curvature of the surface of constant 
t, at the point x. For the metric in our present situation we obtain 

kn = -\e^d t {e 2X ) = -e 2X ^\, 

k 22 = -te-", 

/C33 = —te~^ sin e 0, 

and 

k(t,r) = -e-" + j 
cf. Recalling ( ^ ) we obtain 

2fi\ _ e + e - 2 ^ 



A = e 2 M4vrtp-^-^]>-e 2 ' 1 



y 2t J ~ 2t 
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and 



If e = or e = — 1 then 



and the estimate 



k(t,r)<e^ e 36 ^ . 



k(t,r)<-—e-' t , 



t 

completes the proof of our assertion in these two cases. For e = 1 the estimate 

P-2/S 
g -2 M> ^ >1 = 

~~ t 

which holds for t small, shows that 

W,r)<-—<- w , 

and the proof is complete also in this case. □ 

Since we have already determined the second fundamental form it requires 
little additional effort to investigate the quotients 

1,1 1,2 i,3 
fc l K 2 fc 3 

k ' k ' k 

in the limit t— >0. The question whether these quotients have limits is con- 
nected to the concept of a velocity dominated singularity, cf. Q]. It turns 
out that these limits exist and can be determined in the case of small initial 
data, more precisely: 



Theorem 5.3 Let (/, A,/i) be a solution of (|l.2|) -(|l.6| ) with small initial 
data as described in Theorem 4.1. Then 

t^o k(t,r) 3 

lim MM =lim MM = 2 

fe(t,r) t^o fe(t,r) 3 

uniformly in r G IR. 
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Proof : We have 



k{(t,r) _ tX(t,r) 



k{t,r) t\(t,r) + 2' 
kl{t,r) _ kl(t,r) _ 1 



Prom (^3) we have 



k(t,r) k(t,r) t\(t,r)+2 

p 1 ^ 1 
p 2 2 



As we have seen in the proof of Theorem 4.1 

e 2 M (t,r) < Q t} 

and the estimate on the u;-support of / in that theorem shows that 

p U r )=' / Jl + w 2 + F/t 2 f(t,r,w,F)dFdw<Ct- 3+c 

t z J-p(t)Jo v 

so that 

4vrt 2 e 2/ V(t,r)<C7f ^0, t^O. 

Thus 

t\(t,r)^--, i->0, 

uniformly in r, and the proof is complete. □ 

Remark: As is easily checked a homogeneous and isotropic solution in the 
case e = is given by 

f(t,w,F):=^t 2 w 2 +F), X(t):=\nt, e~ 2 ^ = *ft 2 p(t), 

where cp 6 C 1 ( [0 , oo [) is an arbitrary, nonnegative function with compact sup- 
port. For this solution, the limits considered in Theorem are 1/3,1/3,1/3 
which shows that the result in that theorem need not be true for solutions 
whose initial data violate the smallness assumption. 



6 Going forward in time 

Our main interest in this paper lies the behaviour of solutions as t approaches 
the singularity at t = 0. Nevertheless, it is possible to consider the solutions 
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for £>1. Our main result here is that one can establish a bound on the 
^-support of / and in spite of that the solutions need not exist for all t>l. 
Since this is in sharp contrast to the spherically symmetric, asymptotically 
flat case, cf. Q, and also to well known results for the Vlasov-Poisson and 
Vlasov-Maxwell system we include these considerations. 

If one goes through the proof of Theorem 3.1 , but now for t>l, one 
immediate problem is that (3A) need not define fj, n for all t>l but only 
on some finite time interval. This results in a more restrictive continuation 
criterion in the corresponding local existence theorem the proof of which is 
omitted: 

o o . . 

Theorem 6.1 Let (/, A,/2) be initial data as in Theorem [O]. Then there 



exists a unique, right maximal, regular solution (/,A,/x) of (1.2)~( |1.6| ) with 

o o 

(/, A,/i)(l) = (/, A,/i) on a time interval [1,T[ with Te]l,oo]. If 



and 



supj|u>| | (t,r,w,F) e supp/| < oo 
sup | e 2 ^) | rent, te[i,r[} 



then T = oo. 

The main result of this section is the following theorem: 

Theorem 6.2 Assume that (/, A,/i) is a solution of ([Tp-Q on a right 
maximal interval of existence [1,T[, and 

sup{e 2 ^' r) | reffi,, t€[l,T[}<oo. 

Then T = oo . 

Proof : Assume that T < oo. We show that under the assumption on e 2 ^ we 
obtain the bound 



sup j | iz7 1 | (t,r,w,F) € supp/| < oo, 



which is a contradiction to Theorem 6.1. The proof of the bound on w is 
similar in spirit to the proof of []|, Theorem 2.1]. Define 

P+(t) :=sup{w\(r,w,F) £supp f(t)}, 
P-(t) := m{{w j (r,w,F) £supp/(£)}, 
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assume that P + (t)>0 for some te[l,T[, and let w(t) = w>0 denote the 
w -component of a characteristic in supp/. We have 

w = ^~e 2fl J°° J™ (w^/l + w 2 + F/t 2 -w^/l + w 2 + F/t 2 ^ fdFdw 
l + ee 2fl 

Let us abbreviate 

£ = wyJl + w 2 +F/t 2 -wy/l + w 2 + F/t 2 . 

As long as w(s) > we have the following estimates: If w < then £ < 0. If 
w > then 

_ * 2 (l + w 2 + F/s 2 )-tt; 2 (l + w 2 + F/s 2 ) 
w^fl + w 2 + F / s 2 + w\J l + w 2 + F / s 2 
w 2 {l + F/s 2 )-w 2 {l + F/s 2 ) 



w^l + w 2 + F/s 2 + w\Jl + w 2 + F/s 2 
<C * 



w(s) ' 



and thus 



I f P+(s) rFo C 

(s) < C r— / / wf(s,r,w,F)dFdw + —wU 

sw(s)Jo Jo s 



){s) 

^ C fP + (s) 2 . \ 
s \ w(s) J 

where P + :=max{P + ,0}. Thus 



as long as > 0. Let to G [M[ be defined minimal such that w{s) > for 
se]t ,t]. Then 

/•* 1 

w 2 {t)<w 2 {t ) + C -P + {s) 2 ds. 

J to s 

Now either to > 1 and if (to) = or to = 1 and w(to) < wq. Thus 



P+{t) 2 <wl + C C-P+isfds 
J\ s 
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for all i€[l,T[, since this estimate is trivial if P+(i)<0. If T<oo this 
estimate implies that P+ is bounded on [1,T"[. Estimating w(s) from below 
in case w(s) < along the same lines shows that P_ is bounded as well, and 
the proof is complete. □ 

Remark: In the case of spherical symmetry it is easy to see that a solution 
cannot exist for all t>l, regardless of the size of initial data. For such a 
solution (/,A,/i) the estimate 

t t Jl 

< e-^» + l 
t 

has to hold on the interval of existence [1,T[. Since the right hand side of 
this estimate tends to —1 for t — >oo it follows that T<oo and 

||e 2 ^|Hoo, t->T, 



by Theorem 3.2. This should be compared with the behaviour of the 
Schwarzschild metric in the limit r /* 2M, i. e., approaching the event hori- 
zon from inside of the black hole; cf. our interpretation of the metric ( |1 . 1| ) 
in the introduction. 
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